Controlling high-frequency collective electron dynamics via single-particle complexity 
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We demonstrate, through experiment and theory, enhanced high-frequency current oscillations 
due to magnetically-induced conduction resonances in superlattices. Strong increase in the ac power 
originates from complex single-electron dynamics, characterized by abrupt resonant transitions be- 
tween unbound and localized trajectories, which trigger and shape propagating charge domains. 
Our data demonstrate that external fields can tune the collective behavior of quantum particles by 
imprinting configurable patterns in the single-particle classical phase space. 
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Understanding the interplay between the properties 
of individual objects and their collective behavior is of 
fundamental interest in many fields [1-4 . It explains, 
for example, jamming and pattern formation in granular 
systems (TJ |5]-[7] , dynamical heterogeneity in phase tran- 
sitions [3", tunneling dynamics and quantum phases in 
cold atoms [8, 9 , and the synchronization of networks 
and complex adaptive systems [2l[T0]. Moreover, inter- 
actions between particles play a key role in determin- 
ing the structures formed when the particles come into 
contact [4 . Consequently, tailoring single-particle dy- 
namics may provide a route to controlling the collective 
dynamics of many-body systems. This is a major chal- 
lenge both in fundamental science pTHT3] and for devel- 
oping new technologies such as high-frequency electronic 
devices IHHH], whose performance can be greatly en- 
hanced by applied quantizing magnetic fields [jTl [18] . 

The phase space structure of individual particles, in 
particular the existence and relative location of regu- 
lar and chaotic trajectories, critically affects thermal- 
ization and diffusion both in classical and quantum sys- 
tems [191 120]. Therefore, manipulating the single-particle 
phase space, by generating new chaotic trajectories for 
example, is a promising strategy in the search for ways 
to control other collective phenomena. 

Usually, the transition to chaos in Hamiltonian sys- 
tems occurs by the gradual destruction of stable orbits, in 
accordance with the Kolmogorov-Arnold-Moser (KAM) 
theorem [21 . In far rarer non-KAM chaos, the chaotic 
orbits become abruptly unbounded when the perturba- 
tion frequency attains critical values and map out intri- 
cate "stochastic webs" in phase space [19 . Experimen- 
tal realization of non-KAM classical chaos was recently 
achieved using a quantum system [22 : a semiconductor 
superlattice (SL) with a magnetic field, B, tilted rela- 
tive to an electric field, F, along the SL axis. When the 
frequency of single-electron Bloch oscillations along the 
SL axis is commensurate with that for cyclotron motion 
in the plane of the layers [Ml [22l |23], the orbits map 
out stochastic webs in phase space, which delocalize the 



electrons in real space. This delocalization creates mul- 
tiple resonant peaks in the single-electron drift velocity 
versus F characteristics, which enhance the measured dc 
conductivity [22 . 

In this Letter, we show, via experiments and theo- 
retical modeling, that non-KAM single-particle trajec- 
tories have a dramatic effect on the ac collective dynam- 
ics of electrons in a SL, thus providing a mechanism for 
controlling and enhancing the high-frequency (GHz) ac 
response of this solid-state device. Our work therefore 
overcomes a major long-standing limitation of SL oscil- 
lators, namely how to boost the ac power produced by 
an individual device. At critical applied fields, the single- 
electron phase space reconfigures into extended stochas- 
tic webs, which resonantly increase the electron drift ve- 
locity. This boosts the strength and speed of propagating 
charge domains and, hence, the resulting current oscilla- 
tions. Such tailoring of the collective electron dynamics 
is a non-trivial nonlinear effect, which we explain by re- 
lating islands of high-power current oscillations directly 
to the emergence of web patterns in phase space. 

More generally, our results demonstrate that collective 
effects in periodic quantum systems can be controlled by 
imprinting specific patterns in the underlying classical 
phase space. Such control may be realized in a wide 
range of periodic systems, including quantum cascade 
lasers [TTl [18] and Bose-Einstein condensates in an opti- 
cal lattice [9l[24], or by applying different perturbations, 
for example an ac bias voltage in the case of the SL [25] . 

Our SL was grown by molecular beam epitaxy on a 
(lOO)-oriented n-doped GaAs substrate. It comprises 15 
unit cells [light and dark blue layers in Fig. [TJa) inset], 
which are separated from two heavily n-doped GaAs con- 
tacts [green in Fig. [TJa) inset] by Si-doped GaAs layers of 
width 50 nm and doping density 1 x 10^^ cm~^. Each unit 
cell, of width d and Si doped at 3 x 10^^ cm~^, is formed 
by a 1 nm thick AlAs barrier, a 7 nm wide GaAs quan- 
tum well (QW) and 0.8 In As monolayers at the center of 
each QW. The In As layer facilitates the direct injection 
of electrons into the lowest energy miniband, of width 
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FIG. 1: (Color) (a) dc I{V) curve measured for the SL when 
B = at temperature T = 293 K. Inset: schematic diagram of 
the SL showing the emitter/cohector contacts (green), barrier 
and quantum well layers (dark and light blue respectively), 
and orientation of the applied electric field, F, and magnetic 
field, B, relative to the co-ordinate axes, (b) Color map (scale 
inset right calibrated in dBm) showing the spectrum of /(t) 
oscillations measured over a range of V and frequencies, /, 
near the first harmonic. Left inset: a typical frequency spec- 
trum of the I{t) oscillations measured around the first har- 
monic for y = 0.4 V and T = 293 K. 



A = 19 meV, and also creates a large minigap (A^ = 198 
meV), which prevents interminiband tunneling (see Sup- 
plemental Material [26] and Ref. [27] )• For electrical 
measurements, the SL was processed into circular mesa 
structures of diameter 20 fim with ohmic contacts to the 
substrate and top cap layer. Magnetic field studies were 
performed in a superconductor solenoid magnet with B 
tilted at an angle to the SL (x) axis [Fig. [TJa) inset]. 
An Anritsu Spectrum Analyzer (MS2667C) was used to 
measure the frequency, /, and amplitude of GHz current 
oscillations in the SL. 

Figure [TJa) shows the room temperature (T = 293 K) 
current- volt age, I{V)^ characteristics of the SL measured 
for 5 = 0. The I{V) curve exhibits an ohmic region at 
low followed by a current peak and a region of neg- 
ative differential conductance (NDC) that extends over 
a wide range of applied bias. The NDC creates regions 
of high electron density (propagating charge domains) 
[28] , which generate self-sustained oscillations [291 130] in 
/ versus time t. The inset of Fig. [ijb) shows a typical 
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FIG. 2: (Color) I{V) characteristics (a) measured and (b) 
calculated for the SL at the values indicated when B — 7 T 
and T = 4.2 K. Inset in (b) shows the regions of current insta- 
bility (shaded) in the I{V) curves calculated for ^ = 0° (lower 
trace) and — 40° (upper trace, offset vertically by 20 mA for 
clarity). The maximum (minimum) values of the ac current 
correspond to the upper (lower) edges of the shaded areas. 
Dashed curves in inset indicate the time- averaged current. 



frequency spectrum of the I{t) oscillations, centered on 
the first harmonic and measured for a dc applied voltage 
y = 0.4 V at T = 293 K. The color map in Fig. [ijb) 
(scale inset right) illustrates the evolution of the spec- 
trum as V changes within the NDC region. Red and 
green colors correspond to high signal amplitude. Note 
that the frequency (/i = 1.13 GHz) of the first harmonic 
depends only weakly on V. 

Applying a large (7 T) tilted magnetic field changes 
greatly both the dc and ac current. Figure [2|a) shows 
dc I{y) curves measured for a range of at low tem- 
perature (T = 4.2 K). Increasing from 0° to 60° [blue 
and red curves in Fig. |2|a)] doubles the peak current, 
/p, from ^ 15 mA to ~ 30 mA. Further increasing 9 to 
90° decreases Ip to ^ 25 mA. Figure [2|a) also reveals 
that for ^ 45°, the tilted magnetic field strengthens 
the current instability in the NDC region. This increases 
the amplitude, /^^ of the I{t) oscillations, as shown in 
Fig. ia). 

The inset in Fig. |3ja) shows frequency spectra of the 
I{t) oscillations measured for ^ = with ]/ = 0.3 V (lower 
trace) and 9 = 60° with V = 0.8 V (upper trace) at B 
= 7 T. Note that tilting the magnetic field increases the 
amplitude of the fundamental peaks [marked by arrows 
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FIG. 3: (Color) Color maps showing how the amplitude, 
of the /(t) oscillations varies with V and Q when S = 7 T 
and T — 4.2 K: (a) experimental data - pixellation reflects 
discrete V and values at which we measure /(t) [inset shows 
amplitude of frequency spectra measured for ^ = at V = 0.3 
V (lower trace) and = 60° at = 0.8 V (upper trace, ofl[- 
set vertically by 50 dBm for clarity) when B = 7 T, arrows 
indicate first harmonics]; (b) and (c) are corresponding the- 
oretical plots calculated, respectively, without and with the 
external resonant circuit shown in the inset of (c). 



in Fig. |3|a)] and higher harmonics. To explore this en- 
hancement further, the color map in Fig. [sja) shows 
how the measured la values vary with V and 0. The 
black area indicates (V, 0) values for which / exhibits no 
temporal oscillations. The colored areas reveal two dis- 
tinct islands of large- amplitude (white highest) current 
oscillations, labeled "ET" and "r = 1 & r = 2". For 
30° ^ ^ 75°, Fig. [sj^a) reveals significant enhancement 
of la- The maximum amplitude, attained within the up- 
per large island, is ~ 0.7 mA, which is larger than that 
measured for any V when ^ = 0° 0.3 mA). Note also 
that the range of V over which the SL exhibits I{t) oscil- 
lations [i.e. the width of the colored islands in Fig. |3ja)] 
changes significantly with 0. In particular, this range ini- 
tially decreases as increases from ^ = 0° to ^ ~ 30°, 
but thereafter increases sharply, attaining a maximum at 
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FIG. 4: (Color) Color map showing the single-electron drift 
velocity, Vd, calculated versus F and for B = 7 T. Islands 
of high Vd (whose right-hand edges are regions of negative 
differential velocity) originate either from Bloch oscillations 
(near left-hand dashed line labeled ET) or from resonant cou- 
pling of Bloch and cyclotron motion when r = 1 and 2, which 
occurs along the middle and right-hand dashed curves respec- 
tively. The inset shows a stroboscopic Poincare section of the 
electron trajectories constructed by plotting the crystal mo- 
mentum components (py^Pz) (yellow dots) in the plane of the 
SL layers at integer multiples of the Bloch oscillation period 
Tb = 27V /(jJb when r = 2. Resonances between the cyclotron 
and Bloch oscillations create unbounded web-like structures, 
comprising interconnecting radial and ring-shaped chaotic fil- 
aments as shown in the inset, which are known in the litera- 
ture as "stochastic webs" 19 . 



^ 70°. 

We model the measured effects at T = 4.2 K by making 
self-consistent solutions of the drift-diffusion transport 
equation and Poisson's equation throughout the device - 
see Eqs. (1), (2) and their detailed description in [31]. 
For T < 50 K, temperature changes have little eflFect on 
the electron transport [32 . The solutions give the elec- 
tron density and voltage drop versus t and position x 
throughout the SL (see Supplemental Material, ^6] Fig. 
S3). This model takes into account the field dependent 
electrical conductivity, (7c{B) = (JclW + ('^c'?"c sin ^)^], of 
the contact layers, where (Jc = 3.8 x 10^ l]~^m~^ corre- 
sponds to an electron scattering time, Tc = 90 fs, in the 
contacts. Here, the cyclotron frequency, ujc = eB/m^, 
depends on the effective electron mass, m*, in GaAs. Our 
model also includes the electron scattering time within 
the SL layers, r = ri[Te/{re + t^)]-^/^ = 280 fs determined 
from the elastic interface roughness scattering time Tg = 
38 fs and the inelastic scattering time = 2.1 ps. To 
simulate the measured dc I{V) characteristics, for each 
value of V we time average the calculated I{t) curve. The 
resulting theoretical dc I{V) curves [Fig. [2|b)] are in 
good quantitative agreement with the experimental data 
[Fig. [2](a)] and accurately reproduce both the voltage 
dependence and peak values of / for all 0. 

The inset in Fig. [2|b) shows how the maximum 
and minimum values of the calculated I{t) oscillations, 
marked by the upper and lower edges of the shaded win- 
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dows respectively, vary with V when = 0° (lower curve) 
and 40° (upper curve). In each case, when V is below a 
threshold, Vth^ the maxima and minima of I{t) coincide, 
indicating dc current. But for V > Vth^ these extrema 
separate due to the onset of current oscillations. The 
time averaged current, shown for V > Vth by the dashed 
curves in Fig. [2|b) inset, yields the dc I{V) curves shown 
in the main part of the figure, which agree well with 
the experimental data in Fig. [2|a). Comparing the two 
traces inset in Fig. [2|b) reveals that applying a tilted 
magnetic field significantly increases the amplitude of the 
ac current oscillations. This can be seen more clearly 
from the color map in Fig. (sjb), showing the amplitude 
of the I(t) oscillations calculated versus V and 6>, which is 
similar to the experimental data in Fig. (sja). In partic- 
ular, the model accurately reproduces the enhancement 
of the I{t) oscillations observed for 30° ^ ^ 75° and 
the shapes of the islands in (F, 6) space where oscillations 
occur. 

We now explain the physical reason for the existence 
and relative strength of the two islands of I{t) oscilla- 
tions. At critical values of B and for which the Bloch 
frequency ujb = eFd/h = ruJcCosO, where r is an inte- 
ger, the electron orbits change abruptly from localized to 
unbounded chaotic paths, which propagate through an 
intricate mesh of conduction channels (a stochastic web 
[19]) in phase space [33]: see the inset in Figj4j This delo- 
calization produces resonant enhancement of the electron 
drift velocity, Vd, and, hence, of Ip as increases from 0° 
to 60° . Figure [4] shows a color map of Vd calculated ver- 
sus F and 6 for B = 7 T. The dashed curves show the 
loci of (left to right) the Esaki-Tsu (ET) peak, which oc- 
curs when ujb = and the r = 1 and 2 resonances. 
The ET peak is visible as the lower-left light island in 
Fig. [4| Together, the r = 1 and 2 resonances generate 
a broader region of higher Vd (upper right light area in 
Fig. [4| for between 30° and 75°. The upper and lower 
high Vd islands in Fig. [4] and the negative differential 
velocity (NDV) regions at their right-hand edges, trig- 
ger propagating charge domains, which generate the two 
corresponding islands of enhanced ac power in Fig. [3] 
(a,b). 

The I{t) oscillations are stronger for 30° ^ ^ 75° 
than for 6> = for two reasons. First, Vd is higher near 
the r = 1 and 2 resonances than at the ET peak [Fig. |4] , 
which means that the charge domains move faster and, 
since the NDV is also stronger, contain more electrons 
(see Supplementary Material, [26 Fig. S4). Second, over 
the voltage range corrresponding to the upper island in 
Fig. [4] at different positions in the SL layers F attains 
the values required for both the r = 1 and 2 resonances. 
Consequently, propagating domains associated with both 
resonances coexist in the SL (see Supplementary Mate- 
rial, [26] Fig. S4), thus increasing both the strength and 
higher-harmonic content [as shown in the upper curve 
inset in Fig. |3ja)] of the I{t) oscillations. The striking 



similarity between the islands of high Vd in Fig. |4]and re- 
gions of high-amplitude current oscillations in Fig. [3]^a,b) 
clearly demonstrate that the collective electron dynamics 
are controlled by resonant delocalization of the underly- 
ing single-particle orbits. 

Although the above model explains how the observed 
amplitude of the fundamental I{t) oscillation varies with 

V and 0, it overestimates the measured fundamental fre- 
quency, /i = 1.13 GHz [Figs, [ijb) andjsja) inset], by an 
order of magnitude: specifically, it gives /i ~ 10 GHz. In 
addition, the model predicts that changing V and/or 
can increase the fundamental frequency by a factor ^10 
|3T] , whereas in the present experiment it is almost inde- 
pendent of these parameters. 

To remove this discrepancy, we modified our model to 
include the equivalent reactive circuit of our SL device, 
shown in the inset of Fig. |3jc). Here, C, is a parasitic 
contact capacitance and L, i?, are the inductance and 
resistance respectively of the equivalent circuit. Taking 
C = 2 pF, L = 5 nH [34 , and R = 17 the revised 
model shows that, as in the experiment, the fundamen- 
tal frequency of the current oscillations is ~ 1 GHz for all 

V and 0. Although, for the present sample, the reactive 
circuit produces a ten-fold reduction in the calculated 
frequency of the fundamental peak in the I{t) oscillation 
spectrum, it has much less effect on the amplitude^ la, of 
the oscillations. To illustrate this. Fig. [sj^c) shows how 
la varies with V and when the inset equivalent circuit 
is incorporated in the model. The shapes and positions 
of the islands in this color map are in good quantita- 
tive agreement with experiment [Fig. [sja)], including 
the sudden enhancement of the voltage range and ampli- 
tude of the oscillations when 6 ~ 60°. We attribute the 
smaller amplitude of the measured current oscillations to 
enhanced electron-electron scattering within the charge 
domains. This acts to reduce their electron density and, 
hence, the I{t) oscillations. 

In conclusion, we have shown that the high-frequency 
collective dynamics and I{t) oscillations of electrons in 
SLs can be enhanced by using resonant delocalization of 
single-particle orbits to shape the Vd{F) characteristics. 
We observe this enhancement over a broad range of V 
and values for which the Bloch and cyclotron frequen- 
cies are locally commensurate within the SL. The effects 
of the resonant Vd peaks on the collective charge-domain 
dynamics and resulting I{t) oscillations are insensitive 
to external reactance. Consequently, such effects may be 
exploited to increase the ac power generated by the SL. 
They may also increase the ac frequency if the SL is de- 
signed to minimize its reactance. We emphasize that the 
mechanism for strengthening the current oscillations re- 
lies only on imprinting additional peaks in single-electron 
Vd{F) curves. Since such peaks can be created by apply- 
ing other signals, e.g. ac electrical [25 or acoustic [35] 
waves, the link that we have demonstrated between sin- 
gle and collective transport has broad implication for en- 
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hancing and controlling band transport effects for waves 
in other periodic systems including photonic and acous- 
tic crystals [36l [37], graphene, where relativist ic quan- 
tum scars of unstable periodic classical orbits were re- 
cently identified [38], and ultracold atoms in optical lat- 
tices [i 13121139]. 
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